Is there a symmetry between absorption and amplification in disordered media? 
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Previous studies found a surprising symmetry between absorption and amplification in disordered 
one-dimensional systems once the system size exceeds certain thresholds. We show that this para- 
doxical result is an artifact due to the assumption of a finite output in solving the time-independent 
wave equation, when in fact both the transmission and the reflection amplitude are divergent as a 
result of the large amplification. More sophisticated approaches, such as time-dependent equations 
including field and matter coupling have to be utilized to correctly describe the wave propagation 
of systems with large sizes. 
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Recent observations [y of laser-like emission from dye 
solution emerged with Ti02 nanoparticles have stimu- 
lated intensive theoretical efforts [p|-p^ to investigate the 
properties of disordered media which are optically active. 
From considerations of enhanced optical paths through 
multiple scattering, random systems are expected to pos- 
sess a reduced gain threshold for lasing jl^Jl^. Corre- 
spondingly, one would expect the transmission in disor- 
dered systems to be enhanced with gain. The longer the 
system, the larger the enhancement. Surprisingly, nu- 
merical calculations Q based on time-independent wave 
equations showed that for large systems the wave propa- 
gation is still attenuated, indicating a localization [ [l5|Jl6| 
of waves even with gain. The rate of the exponential 
attenuation is the same as if the system was absorb- 
ing. Such a symmetry between absorption and amplifi- 
cation was subsequently shown jl^] to hold for the time- 
independent wave equation. 

Intuitively one would expect that the presence of am- 
plification should facilitate wave propagation, not sup- 
press it, even in disordered systems. The peculiar result 
of reduced transmission in gain media is also absent when 
wave propagation in disordered media is described by 
time dependent diffusion equations which always predict 
an increased output and a gain threshold above which the 
system become unstable []l3| , |l4| , p^ . The diffusion equa- 
tion neglects the phase coherence of the wave and is ad- 
equate only when the wavelength is much smaller than 
the mean free path. Thus it was proposed that the 
paradoxical phenomenon may indicate enhanced localiza- 
tion due to the amplification of coherent backscattering. 
However, amplification of backscattering does not neces- 
sarily imply a reduction in transmission since no conser- 
vation of the total photon flux is required in gain media. 
To fully understand the origin of this non-intuitive re- 
sult, we will examine the validity of the solutions derived 
from the time- independent wave equation which has been 
commonly employed in describing the wave propagation 
in active media. 

Linearized time-independent wave equations with a 
complex dielectric constant have been successfully uti- 



lized to find lasing modes by locating the poles p9| ] 
in the complex frequency plane and to investigate the 
spontaneous emission noise below the lasin g th reshold in 
distributed feedback semiconductor laser s |po| , pT|| . Such 
equations are known to be inadequate [ p9| , p0[ to describe 
the actual lasing phenomena due to their simplicity in 
dealing with the interactions between radiation and mat- 
ter. However, it is generally believed that the time- 
independent equation should suffice to describe the wave 
propagation in amplifying media, before any oscillations 
occur. We will unambiguously show that the so-called 
symmetry p7[ | between amplification and absorption is 
an artifact due to the unphysical assumption of a fi- 
nite output in solving the wave equations. We show 
that for each system, there is a frequency-dependent gain 
threshold above which both the total transmission and 
the total refiection become divergent. Solving the time- 
independent wave equations by incorrectly assuming a 
fixed output leads to unphysical solutions that does not 
correspond to the true behavior of the system. 

To demonstrate our point, we take the simplest exam- 
ple of a uniform active media sandwiched between two 
mirrors as feedback (see the insert in Fig. 1), the clas- 
sical Febry-Perot setup. Wave propagation within the 
active media is simply described with the following phe- 
nomenological wave equation, 

where E{z) is the electric field and the dielectric constant 
e{z) — e'{z) — ie (z) with the imaginary part signifying 
amplification (e > 0) or absorption (e < ). We 
point out that in electromagnetic theory, the imaginary 
part of the dielectric constant is proportional to the con- 
ductivity of the material and thus cannot be negative. 
The negative dielectric constant is strictly speaking only 
an effective way to introduce coherent amplification [g2| . 
Complex potentials known as optical potential have long 
been employed in nuclear physics to describe nuclei scat- 
tering processes. 
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The transmission and the reflection amphtude can be 
easily obtained by solving Eq.(l) to yield: 



where ti — 2fc/(fc + feo), t2 = 2fco/(fc + fco), and ri = r2 — 
{k — ko) / {k + ko) are the transmission and reflection coef- 
ficients at the left and right two interfaces, respectively. 
ko = V^ox and k = k' - ik" = y/e' - ie"^ are the 
wavevectors outside and inside the system, L is the dis- 
tance between the mirrors (system size). The oscillation 
condition for lasing is correctly given by 1 — rir2e^*'^^ = 
at which both the transmission and reflection coefficient 
diverge. 

However, Eq.(2) also predicts the exponential decrease 
of the transmission coefficient for large system sizes. In 
fact the term in the numerator with the exp{ikL) in- 
creases exponentially as the length of the system L in- 
creases because of the gain, but the term in the denom- 
inator with the exp(2ikL) increases even faster, mak- 
ing the transmission coefficient decay asymptotically as 
\t\ ~ \exp{—ikL)\ = exp{—k"L). This is clearly shown in 
Fig.l, where we plot the ln{T) versus L for a one-layer 
system. Notice that for large L, ln{T) decreases as L in- 
creases despite the fact that we have gain at here. Gain 
effectively becomes loss at large lengths! Remember the 
system is homogeneous thus disorder is definitely not re- 
sponsible for this strange behavior. Thus the inhibition 
of wave propagation for large systems is clearly not a 
result of amplification in backscattering. 

A clear picture of what is going on can be obtained 
from the path integral method In such an ap- 

proach, the total transmission coefficient can be obtained 
by adding the paths of the successive reflected and trans- 
mitted rays. In doing so, we obtain that 

t = iit2e^'=^[l + (rirae^''^^) + {r^r^e^^'^'^f + ■ ■ ■] (3) 

where the first term represents the direct transmission of 
the incoming wave, and the second term represents the 
wave which was reflected first by r2 at right interface and 
then by ri at left interface and subsequently transmitted 
through. More terms from sequences of multiple trans- 
missions and reflections from the two mirrors follow. It 
is clear Eq.(3) represents an infinite series whose sum re- 
produces Eq. (2), provided that the following condition 
is met, 

|rir2e2''='^e2'="^| < 1. (4) 

When the condition given by Eq.lQ) is violated, such 
as when the system size or the gain is large, the physical 
output represented by the sum diverges, even away from 
the oscillation pole. 

We hope to note at here that the sum of right side of 
Eq.(3) include the phases of all paths, so that it includes 
all interference of scattering waves. We also note that the 



results of Eq.(3) are consistent with time-dependent the- 
ory because different terms of Eq. (3) can explain as out- 
puts at different time from a same incidence or the output 
at same time from a series incidence at different time. By 
this way, the out-put of a system with gain much larger 
than the threshold will increase exponentially versus time 
after a plane wave incidence, but time-independent the- 
ory suppose a small output as if the system is of absorbing 
material. Base on time-dependent Maxwell equations, 
a well-developed FDTD (finite-difference time-domain) 
method can help us to see the obvious different behav- 
ior of transmission of two kind of systems versus time. 
We choose a plane wave with A = 800nm as the inci- 
dence from left to the setup of Fig. 1 with L = 4300nm. 
The dielectric constant is taken to be £o = 1 outside and 
Re{e) = e' = 9 inside. The threshold of this system is 
e" = 0.12 from Eq.(4). Our numerical results show that, 
after a short instantaneous state, the output of system 
will get to a stable value which is the transmission of 
the system if the system with under-threshold gain or 
with absorption. Unless the gain or the length is over 
the threshold, the results of time-independent theory are 
same as the time-dependent theory. But if the system 
with gain is above the threshold, the output of system 
will increase exponentially versus time as predicted by 
Eq.(3). In Fig. 2, we plot the logarithm of amplitudes of 
out-put at right side of system versus time for different 
gain, we can see the numerical results exactly same as 
we predicted. We also examined the slope of lines with 
the gain larger the threshold in Fig. 2, and found the 
slope is exactly same as we predict in Eq. (3), equal to 
logw\rir2e'^''" ^\. 

The divergence of the transmission above threshold 
even away from the oscillation pole is the key in under- 
standing the failure of Eq.(2), which is conventionally 
derived from boundary condition by implicitly assuming 
that the output is finite. Normally the physical bound- 
ary condition is satisfied when tii2e**'^ -f rir2e^*'^^t = t, 
resulting in Eq. (2). Obviously this condition lost its 
meaning when t is infinite. As a result, Eq. (2) becomes 
unphysical. For a lossy system, the condition given by 
Eq.(^ is always satisfied. Thus the divergence is a new 
phenomenon occurring only to systems with gain. The 
breakdown of the time-independent wave equation sig- 
nals the large fluctuations of the transmission in time 
and calls for more sophisticated theories that can take 
into account the interaction between radiation and mat- 
ter to correctly describe the response of the system. 

To show that the exponentially attenuated results 
of time-independent theory with gain is not from the 
backscattering effect, but from a theoretical mistake, we 
also do the same numerical experiment to random sys- 
tem with gain. The system is made of 50 cells, each cell 
contains two kind layers with dielectric constant ei = 1 
and £2 = 4 — ie" . To introduce disorder, we choose the 
width of first layer of the nth cell to be random variable 
a-n — ao(l + W^7),where oq — 95 nm, W describes the 
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strength of randomness which is 0.8 in our model and 
7 is a random number between (—0.5,0.5). The width 
of second layer in nth cell is 6„ — 215nm — a„. The 
wavelength of incident wave is A = 1200nm. In Fig. 3, 
we plot the logarithm of output at right side of system 
versus time with plane wave incidence for different gain. 
We also get the same results as what we predicted, the 
exponential increase of the output versus time when the 
gain is over the threshold value, which are showed in 
Fig. 3. The time-dependent results also show that the 
threshold of random system is quite small as previously 
predicted and showed experimentally [||. Actually all 
time-dependent theories about the system with the gain 
or the length over the threshold, such as Letokohov dif- 
fusion theory Lamb theory give the divergent 
output of the system. The physical meaning of the di- 
vergent output is that the rate of generating photons by 
the induced radiation is larger than the escaping rate of 
photons from the interfaces of the system. Now the field 
in the system will become stronger and stronger even if 
the frequency of the wave is not resonant frequency(at 
the resonant frequency, the system become a laser). 

For multilayer systems, to see whether the results 
of transmission or reflection coefficients of the time- 
independent theory is physical, a simple approach is to 
check every layer of the system by the following method. 
When we examine one certain layer, such as m-th layer, in 
the multilayer system, we can assume that the left part 
of system forms an effective interface of the layer with 
transmission and reflection coefficients ti and r; , and the 
right part forms the other effective interface with tr and 
rr- Here we assume both right and left subsystems are be- 
low threshold. The convergence condition for every layer 
of the system is then given by Eq.(2), with (ti, ri), and 
(^2j ^2) replaced by (ti,ri) and {tr, r^), respectively. For 
every layer, Eq.(2) gives us a line. The cusps are formed 
by the lowest limit out of all lines. Fig. 4 shows the re- 
sults of such a calculation for a 40 layer system, which 
is formed by two kind of layers with widths Li — 95nm, 
L2 = 120nm and dielectric constants e'l = 1, £2 = 2 re- 
spectively . In reality, the curve only gives an indication 
of the magnitude of the gain above which the results from 
time-independent equation becomes suspect. 

In solving the time-independent equations, it is dif- 
ficult to know exactly when the solutions break down. 
One certainly cannot tell from the expression of the total 
transmission and the total reflection coefficient which are 
well behaved except exactly at the oscillation pole. How- 
ever, a rough indication is that when the gain for a fixed 
system approaches the lasing threshold for nearby poles 
or when the system sizes exceed the threshold length, the 
solutions should not be trusted. However, the threshold 
condition is still given by the poles in the complex plane. 
A study of the distribution of these poles has been car- 
ried out through carefully locating the pole position by 
continuously tuning the gain up at a fixed frequency until 
the transmission coefficient diverges |8| . 

The implication of the analysis above is that the cal- 



culation of transmission and reflection coefficient with 
the traditional method become suspect once the gain 
or system size reaches a certain value. Care has to be 
taken to ensure that the system is not above threshold 
and the solution is physical. This remark, unfortunately, 
also pertains to the application of the powerful invariant 
embedding method |P,|7|,p|,p^ . in which the finite reflec- 
tion and transmission are assumed implicitly. 

In summary, this note aims to bring attention to some 
peculiar aspects of the time-independent wave equation 
when the gain is above certain threshold for a flxed 
length, or equivalently, when the length of the system 
exceeds a certain value for a fixed gain. Thus some of 
the conclusions on the statistical properties of the re- 
fiection and the transmission coefficient in media with 
gain become suspicious at large system sizes. The time- 
independent equation is inadequate to describe the am- 
plification of light under these conditions. Nevertheless, 
the simplicity of the time-independent equation can be 
used effectively to locate resonant conditions, even in dis- 
ordered systems. A complete treatment of the wave prop- 
agation in gain media may require the construction of 
the time-dependent solution out of the continuous and 
discrete solutions of the time- independent equation [ p5| . 
Unlike for the Hermitian system, the completeness of 
these solutions could not be proved easily ||26| when the 
potential is complex. 
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FIG. 1. The logarithm of the transmission coeffi- 
cient, ln{T), versus the length L of a simple Febry-Perot 
setup. The dielectric constant of the Febry-Perot device is 
e = e' — is" = 2 — iO.Ol. The dielectric constant of outside 
medium is e = 1 without gain. 

In the insert a schematic representation for wave transmission 
in a simple Febry-Perot setup is shown. 

FIG. 2. The logarithm of output at right side of system 
versus time of a Febry-Perot setup with a plane wave (A = 800 
nm)incidence. The length of system is the L = 4300 nm and 
the dielectric constants of inside and outside materials are 
£i = 9 — ie" and £1 = 1 respectively. The time unit to is 
the time for the wave to travel back and forth through the 
system. The critical gain of the system is e" = 0.12. 

FIG. 3. The logarithm of output at right side of system 
versus time of a random system with a plane wave(A — 1200 
nm) incidence. The system is made of 50 cell and each cell 
contain two kind layers with dielectric constant si = 1 and 
£2=4 — ie" respectively. In nth cell, the width of first layer 
is a„ = a()(l + Wj), whore ao = 95 nm, W = 0.8 and 7 is 
a random number between (—0.5,0.5), the width of second 
layer is bn = 215 — a„ nm. The time unit to is the time for 
the wave to travel back and forth through the system with 
phase speed. 

FIG. 4. The negative value of the imaginary part of the 

dielectric constant e" versus the incident wavelength A for a 
40-laycr system with two kind of layers of width Li = 95nm 
and L2 = 120nm, and with dielectric constant s'l = 1 and 
s'l = 2 respectively. The solid line shows the maximum value 
of e" which the solution of the time-independent equations 
can be trusted. 
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